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Abstract 

An ab initio theory for Fano resonances in plasmonic nanostructures and metamaterials is de- 
veloped using Feshbach formalism. It reveals the role played by the electromagnetic modes and 
material losses in the system, and enables the engineering of Fano resonances in arbitrary geome- 
tries. A general formula for the asymmetric resonance in a non-conservative system is derived. 
The influence of the electromagnetic interactions on the resonance line shape is discussed and it is 
shown that intrinsic losses drive the resonance contrast, while its width is mostly determined by 
the coupling strength between the non-radiative mode and the continuum. The analytical model 
is in perfect agreement with numerical simulations. 
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1. INTRODUCTION 



Asymmetric resonances display unique features, compared to their symmetric Lorentzian 
counterpart, and are currently the subject of considerable research efforts in photonic and 
plasmonic nanostructures. A theoretical derivation was first proposed by Fano to explain 
autoionization of atoms and the asymmetric shape of these resonances that now bear his 
name 1 . In fact, the interference phenomenon underlying Fano resonances is a general wave 
phenomenon, appearing in particular as Wood anomalies in gratings^, in extraordinary op- 
tical transmission 3 , dielectric 4 and metallic 5 photonic crystals, and more recently in op- 
tomechanical systems^, plasmonic nanostructures^ - —, or as the plasmonic analog of electro- 
magnetically induced transparency (EIT) in metamaterials 1 ^— . Fano resonances exhibit 
a very strong sensitivity to changes of the local environment as well as a sharp spectral 
dispersion. The publications of Luk'yanchuk et al.— and Miroshnichenko et al— provide 
extensive reviews on the state of the art of Fano resonances in nanoscale structures and 
their future application prospects, which include sensing, optical modulation and switching 
or non-linear devices. Although Fano resonances have been studied in a broad variety of 
complex plasmonic nanostructures and metamaterials, their analysis relies either on a clas- 
sical oscillator modet^* 1 ^, phenomenological models 4 -* 1 ^, coupled-mode formalism^ or the 
quantum mechanical approach used by Fano 1 - to fit experimental data and to understand 
the mechanisms behind the resonance shape. For example, it is not well understood in 
a realistic plasmonic system composed of several interacting particles how the individual 
modes and their coupling affect the overall Fano-like resonance of the system. The fact that 
a discrete radiative mode supported by a plasmonic nanoparticle can act as a continuum in 
Fano interferences still remains a challenging theoretical task. 

We develop in this work an ab initio theory for asymmetric resonances for electromagnetic 
scattering in general dispersive and lossy media. A general formula for the asymmetric 
resonance in a non-conservative medium is derived [Eq. ( I24p and Fig. [1] . This theory reveals 
the role played by the electromagnetic modes and material losses and enables the engineering 
of Fano resonances. The influence of electromagnetic coupling onto the resonance line shape 
is illustrated with the numerical example of a realistic plasmonic nanostructure. It is finally 
shown that our final result, Eq. (j24p . is also valid for the mechanical model of two coupled 
oscillators. 
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FIG. 1: Resonance shape function a /a as a function of the reduced frequency k = (u 2 — — 
a^A) /r for different values of the asymmetry parameter q and the screening parameter b [Eq. ()24|) ] . 
Intrinsic losses affect the Fano resonances contrast described by the parameter b, induce a spectral 
broadening and a damping of the asymmetry parameter. 

2. ELECTROMAGNETIC THEORY OF FANO RESONANCES 

We consider the electromagnetic scattering on a dielectric or metallic object in a dielectric 
background. The scatterer supports a non-radiative (dark) mode interacting with a contin- 
uum or a radiative (bright) mode. A harmonic time-dependence for the fields E = Eoe~ lu)t is 
assumed throughout. The relative dielectric permittivity e(r, uj) is generally complex and has 
frequency dispersion. The frequency dependant electric field E satisfies the wave equation 



To simplify notation, in the following the dependance of the electric field E on u is kept 
implicit. Let us assume that the scatterer permittivity is given by a Drude model with 
plasma frequency u p ; the scaling law of Maxwell's equations allows then to scale all frequency 
units by u p and length units by 2ttc/u p . We introduce the generally complex and frequency 
dependant differential operator M defined by 




(1) 




(2) 



The wave equation can be written for a vectorial wave function |E) 



CM -w 2 I)|E) =0, 



(3) 



where I is the identity operator. The inner product is defined by 




(4) 
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Following Feshbach, Bhatia and Temkin, we introduce the orthogonal projection operators 
P and Q splitting the wave function into a radiative (bright) and a non-radiative (dark) 
partj22>2i. Any wave function |E) can be decomposed as |E) = Q\E) + P|E) where only P|E) 
satisfies the radiation condition. Equation becomes 

(M^-co 2 l)(Q\E)+P\E}) = 0, (5) 

yielding the two coupled equations: 

(QM^Q ~ uj 2 1)Q\E) = -QM^P\E) (6) 
(PM^P - u 2 l)P\E) = -PM^Q\E) . (7) 

We consider a unique non-radiative mode |Ed), defined to be eigenfunction of the projector to 
non-radiative modes Q\E d ) = \E d ), and to satisfy QM t Q\E d ) = zj\E d ) and |(E rf |E d )| 2 = 1. 
Taking material losses into account, the quantity z d = u d + i^ d is generally complex. Its 
real part u d is the mode resonance frequency, and ^ d its intrinsic damping. The system is 
studied in the vicinity of the resonance frequency u d . In this frequency range, we assume 
that the space of eigenfunctions of Q is spanned by \E d ) only^ 2 -, so that 

Q=\E d )(E d \. (8) 

Equation (J6]) can be written in Q-space as: 

Q\E) = —L^\E d )(E d \M P|E) . (9) 
u 2 - zj — <" 

Inserting Eq. into Eq. ((7]), a wave equation with source is obtained for P|E): 

(PM P - u 2 l)P\E) = 1 PM \E d ){E d \M P\E) . (10) 

The radiative (bright) wavefunction |PE&) is defined to satisfy Eq. (fTOl) if the mode \E d ) is 
removed. The interference between |PE&) and \E d ) will result in a Fano-like resonance for 
the total wavefunction |E). From its definition, |PEb) satisfies the following homogeneous 
equation: 

(PM^P - u 2 )\PE b ) = . (11) 

Equation (|T0l) can be solved using the dyadic Green's function G_ b of Eq. ( II ip : 

{E d \M PIE) 

P|E) = P|E 6 ) + ^ 2 GPM\E d ) . (12) 
z d — cu 2 — b — w 
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In Eq. (fl2]l . the wavefunction |E) is related to |E) but does not have the same asymptotic 
behavior; an expression for |E) is derived later. Multiplying Eq. ( 1T2|) on the left by (Erf|Aj, 
and replacing (E d \M P|E) in Eq. §F2§ yields: 

(E d \M \PE b ) 

PW = l PV i) + ^—^ g b PM >{ V d) , (13) 

where 

A = -(E d \M^Pg h PMjE d )/u d (14) 

is the shift in the resonance position u d due to the field overlap between the continuum 
|PE^) and |E^). Using Eqs. (|9]) and ( |T2|) . we have a final expression for |E): 

(E d \M \PE b ) 

IE) = |PE b ) + \, ' ' (g^E,) " IE,)) • (15) 

The wavefunction |E) must have the same norm as P|E{,) in the far-field (i.e. |P|E)| 2 = 
|P|E&)| 2 ). The asymptotic behavior of |E) is now compared to \PE b ), leading to a relation 
between |E) and |E). Most of plasmonic nanostructures and metamaterials are embedded in 
a dielectric medium, either homogeneous or on a substrate, allowing us to assume that the 
permittivity is real and non dispersive in the radiative region. Therefore P.M is self-adjoint 
and the set of solutions |PE{>) of Eq. (|TT|) forms an orthogonal basis of modes. The Green's 
function G_ b is expanded on this continuum : 



yielding: 

\(E d \M \PE b 



P|E) = |PE fe ) 



2tu(z 2 -u 2 + cu d A) ' 
Let us define the intrinsic damping parameter: 



(17) 



l(E rf |^|PE,)| 2 7d ^ 
1 Lo(co 2 -co 2 + co d A) 2 ' [ } 



the resonance width: 



(E d |^|PE b 
2w(l - Ti) 



(19) 
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and the reduced frequency: 

k= (u 2 -Lu 2 d -co d A)/T. (20) 

In order for |E) to have the same normalization in the far-field as \E b ), it is related to |E) 
by 

|E) = cos7 ? |E)/(l-r i ), (21) 

where cotanr/ = k, considering 7^ <C u d and neglecting any second order contribution. We 
have asymptotically P\E) = P\E b ) exp(ir]), meaning that the phase of |E) shifts rapidly by 
~ 7r in a frequency region T around the resonance position u d + u d A, resulting in a Fano- 
like interference. The parameters u d and A determine its spectral position. Its width T is 
influenced by the field overlap between the continuum \PE b ) and the non-radiative mode 
I Erf), which is a direct effect of their interference. We consider now a transition operator T 
from an initial excited state \g) to a final state |E) with emission of a photon, which can 
be interpreted as the response of an optical system to an external excitation, such as the 
Local Density of States (LDOS)^, the forward scattering cross-section^ or the reflectance 
of a two-dimensional array. From Eq. f[T5]) : 

( 9 \T\E) \{E d \M^\PB b )\ 2 

1 O, .(Jl . .2 , , . AW1 \1 ' 



(g\T\PE b ) 2u(zj -u* + u d A)(l - r { ) 

with 

(g\T\E d ) - (g\Tg b PM \E d ) 
g = 2w(1 " ri) ((E d \M^\PE b )y(g\T^E b ) ■ (23) 

The parameter q is given by the ratio between the optical response of the perturbed non- 
radiative mode and the continuum. Assuming q is real, one gets from Eqs. (j2ip and ( 122]) 
the ratio of the optical response of the total field |E) to the one of the continuum \PE b ): 

\(g\T\E)\ 2 U + gf + b 
* ~ \(g\T\PE b )\* ~ a ^ + 1 ' [M) 

where b = K 2 r^g 2 (l — T^) -2 and a = (1 — r.;)~ 2 . Equation ff24|) is a generalization of Fano 
formula 1,21 to vectorial fields and lossy materials. It describes in particular the resonance 
strength in plasmonic nanostructures. If the non-radiative mode has no intrinsic losses, one 
has Ti = 0, b = and the Fano formula is recovered. The parameters T, Tj, q, A and b in 



6 



Eq. (|24|) are assumed to be constant in the resonance region, but one could consider their 
lowest energy- dependant corrections over a larger range of frequencies^. 

The parameter b does not play a role for the spectral width of the resonance but it 
influences its contrast (Fig. [T]), preventing it to reach zero values. Assuming the parameter 
b is constant over a frequency region around ojd, it can be evaluated at u = u^, leading to: 



counting for the overlap between the non-radiative mode |E^) and the continuum |PE&). If 
the coupling between the two modes is too weak compared to intrinsic losses, the parameter 
b increases, which has for effect to screen the Fano resonance (Fig. [I]). The screening pa- 
rameter b is a limitation for the resonance width T, which appears to be the main effect of 
intrinsic losses in Fano resonances. It is also a critical parameter for achieving electromag- 
netically induced transparency since it determines the maximal amount of light that can be 
transmitted^ 1 ^. Intrinsic losses included in Tj are also responsible for a damping of the res- 
onance width r and the asymmetry parameter q. The parameters q and b together describe 
the resonance shape, and depend on what the transition element T is chosen to represent: 
for instance the local density of states (LDOS) 23 , the forward scattering cross-section^ or 
the reflectance of a two-dimensional array. 

3. NUMERICAL VALIDATION 

We now verify numerically the validity of Eq. (1241) with the example of a dolmen-type 
plasmonic structure^^, made of three metallic beams arranged as in Fig. |2j The method 
used to model electromagnetic scattering on plasmonic structures is based on surface inte- 
grals^ - —. In the quasistatic approximation, a dark mode cannot be excited by a planewave, 
but by a rapidly varying field, such as a dipole placed in the near-field of the structure^. 
The local density of states (LDOS), defined as the imaginary part of the Green's tensor 
at the position of the dipole, provides all the spectral and scattering information for the 
structure^. It is computed in Fig. [2] at one end of the double beam structure, revealing the 
existence of a non-radiative quadrupole mode at an energy of hwd = 1.14eV. The dipolar 
mode of the perpendicular beam has a larger spectral width due to radiative losses and 



b ~ 4 




(25) 



In the denominator of b appears the product 
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FIG. 2: Structural decomposition of a Fano resonance in a dolmen-type plasmonic structure in 
air: beam 1 supports a radiative dipolar mode, whereas beams 2 and 3 support a non-radiative 
quadrupolar mode, (a) Local density of states (LDOS) of a dipole emitter placed 50nm from the 
end of one beam of the double beam structure, (b) reflectance of a single beam, (c) reflectance 
of the composite system. Insets: meshed structures and geometry parameters; Iq =160 nm, 1% = 
I2 =300 nm, w =80nm and t =40nm and g = 30nm. For the calculation of reflectance spectra, 
the structures are arranged on a two-dimensional array with period 500nm to avoid neighbor-to- 
neighbor interactions. A Drude model with plasma frequency lo p = 1.367 x 10 16 s _1 and damping 
oj-y = 0.009wj, has been chosen for the scatterers' material. 



ensures the coupling of the non-radiative mode to the radiative continuum. The reflectance 
spectrum of the dolmens array is calculated and divided by the reflectance of an array of 
single beam in order to obtain the shape function a shown in Fig. [TJ In Fig. [31 the distance 
g between beam 1 and the pair of beams 2 and 3 is varied, which changes the strength 
of the overlap and the coupling between the non-radiative mode and the continuum. For 
such a realistic plasmonic system for which no analytical solution exist, its optical response 
does satisfy Eq. (|24p . This equation is in perfect agreement with numerical calculations, 
and allows us to study and engineer its different parameters. As expected from Eq. ([TBI , 
the resonance shift A decreases with decreasing coupling , along with its width T, and the 
absolute value of the asymmetry parameter q decreases. The coupling strength is therefore 
a critical parameter for engineering the width of a Fano resonance. In Fig. HI the Drude 
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FIG. 3: Influence of the gap width g on the different parameters of the resonance for the dolmen- 
type structure in Fig. [2j simulation and fit to Eq. (|24p . Each point in the bottom panels corresponds 
to the fitting parameters extracted from a full calculation for a specific gap distance g. 
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FIG. 4: Influence of the Drude damping cu 7 on the different parameters of the resonance for the 
dolmen- type structure in Fig. [2j simulation and fit to Eq. (|24p . Each point in the bottom panels 
corresponds to the fitting parameters extracted from a full calculation for a specific Drude damping 
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FIG. 5: (a) Mechanical model of Fano resonances: two coupled oscillators of resonance frequencies 
ujb and ujd, and damping 75 and jj- One of them is forced by an external excitation of amplitude /. 
The coupling constant is g. (b) Amplitude of the forced oscillator as a function of the excitation 
frequency for u>b = 50., uid = 50.1, g = 2.0, 7& = 0.1 and 7^ = 0.01. 

damping o; 7 of the beams 2 and 3 is varied, affecting the imaginary part of the non-radiative 
mode's eigenvalue 7^ Of particular importance is the invariance of the resonance shift A, 
in agreement with the theory which has shown that it is independent on 7^. The resonance 
width T increases with increasing losses, along with the b parameter; and the q parameter is 
damped. Overall, intrinsic losses affect the resonance contrast. In order to engineer a sharp 
and contrasted Fano resonance, an optimal value for both V and b has to be found. 



4. MECHANICAL ANALOG: TWO COUPLED OSCILLATORS MODEL 

In this Section, we show that the final result Eq. ( |24|) . that we have derived from electro- 
magnetic considerations, can also be obtained for the case of two harmonic oscillators with 
losses. Consider the two harmonic oscillators system in Fig. [5] which model classically the 
Fano resonance^ 1 ^. The first oscillator modeling the bright (radiative) mode has a reso- 
nance frequency u;& and high damping 75 representing radiative losses. The second oscillator 
modeling the non-radiative (dark) mode has a resonance frequency Ud and low damping 7^ 
The two oscillators have a coupling g. The highly damped oscillator is driven by a harmonic 
external force with amplitude fe lult , representing the coupling to an external field. The 
equations of motion can be written as: 

%b + Ibib + Ub x b + gxd = fe tuJt (26) 
Xd + laid + uj^d + gx b = . (27) 
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The displacement Xj, and Xd of the oscillators is harmonic, therefore x& = Cbe tujt and Xd = 
Cde tujt . The amplitude of the highly damped oscillator is given by: 

r _ <4 ± lldUJ ~ ^ f (9R\ 

b / 2 i * 9\/ 2 i ■ 9\ 9J 1 l z,c V 

(w£ + t'jbuj - uj 2 ){ujd + l ld,u - tJ 2 ) - g 2 

We now develop the denominator of Eq. f l28|) around Co^, the frequency around which the two 
oscillators interfere, giving rise to an asymmetric resonance. Considering that 7d 7& 
Ub, ujd, the quantity u 2 + i^b^ — u 2 is slowly varying in a small frequency interval around Ud 
and can therefore be considered as constant [Fig. E^b)] . It is evaluated as C = u 2 +i r ybOJd—uj^. 
Around Ud, Equation becomes: 

o^ rl ti +i ^~i J - (29) 
C(uj + 27^0; - w 2 ) - g 2 

Neglecting second order contributions from 7^, one obtains from Eq. ( |29i) : 

|c fe | 2 _ (K + q) 2 + b 

W~ k 2 + 1 ' w 

where |A| 2 = |/| 2 /|C| 2 is the amplitude of the forced oscillator if there is no coupling; 
k = (u 2 — u^ — ujdA)/T is the reduced frequency; A = [{uj^ — u 2 )g 2 }/ (\C\ 2 Ud) is the resonance 
shift; T = •jbUdg 2 /\C\ 2 its width; q = (a;J — u>^) / (jb<^d) the asymmetry parameter and 
b = 7j|C| 4 /(7fefi' 4 ) the screening parameter. The latter also satisfies b = 7jg 2 /A 2 which is a 
similar relation between the resonance parameters as in Eq. (|25|) . Hence, Equation (|30|) is 
the mechanical equivalent to Eq. ff24]) . 

The mechanical and the electromagnetic models both feature the same four independent 
parameters A, T, q and b describing the asymmetric resonance. The parameter b appears 
if the intrinsic damping 7d 7^ and increases if 7^ increases, having for effect to screen the 
resonance (Fig. [I]), i.e. decrease its contrast. It should also be noted that an increase of 
the coupling g leads to an increase of the resonance width T, an increase of the absolute 
value of the shift A and a decrease of the screening parameter b, which is in agreement with 
the electromagnetic theory. The parameters A and q are function of the relative spectral 
positions of the two oscillators' resonances a; 2 - — u 2 and vanish if ojf, = Ud- 

Finally, let us mention that this section shows that Eq. f )24l) is not only valid in electro- 
magnetism, but also for any physical system that can be modeled by two classical damped 
coupled oscillators. It is therefore considered as a general formula for the asymmetric reso- 
nance in a non-conservative system. Unlike the electromagnetic theory, Equation f )30|) can 
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not provide an expression of the shape parameters q and b including the detection method 
(for instance back-scattering or forward cross section). This observation appears as a limi- 
tation of the mechanical model, compared to the electromagnetic theory. 

5. SUMMARY 

We have derived an ab initio theory for Fano resonances in plasmonic nanostructures 
and metamaterials using Feshbach formalism. The influence of the electromagnetic inter- 
actions on the resonance line shapes has been discussed. The novel formula that describes 
the resonance line shapes introduces a critical screening parameter driven by the intrinsic 
losses. The resonance width is mostly determined by the coupling strength between the non- 
radiative mode and the continuum. The formula can be derived from the classical model of 
two damped coupled oscillators and is in perfect agreement with numerical simulations of 
complex plasmonic systems. A deep insight into the study of Fano resonances in metallic 
nanostructures has been enabled. 

Funding from CSEM and CCMX-Fanosense as well as stimulating discussion with M. 
Schnieper and A. Stuck are gratefully acknowledged. 
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